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Abstract: We extend the previously obtained results for the thermodynamic potential of 
hot QCD in the limit of large number of fermions to non- vanishing chemical potential. We 
give exact results for the thermal pressure in the entire range of temperature and chemical 
potential for which the presence of a Landau pole is negligible numerically. In addition 
we compute linear and non-linear quark susceptibilities at zero chemical potential, and the 
entropy at small temperatures. We compare with the available perturbative results and 
determine their range of applicability. Our numerical accuracy is sufficiently high to check 
and verify existing results, including the recent perturbative results by Vuorinen on quark 
number susceptibilities and the older results by Freedman and McLerran on the pressure 
at zero temperature and high chemical potential. We also obtain a number of perturbative 
coefficients at sixth order in the coupling that have not yet been calculated analytically. In 
the case of both non-zero temperature and non-zero chemical potential, we investigate the 
range of validity of a scaling behaviour noticed recently in lattice calculations by Fodor, 
Katz, and Szabo at moderately large chemical potential and find that it breaks down rather 
abruptly at fi q > ttT, which points to a presumably generic obstruction for extrapolating 
data from small to large chemical potential. At sufficiently small temperatures T <C fJ, q , 
we find dominating non-Fermi-liquid contributions to the interaction part of the entropy, 
which exhibits strong nonlinearity in the temperature and an excess over the free-theory 
value. 
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1. Introduction 

At large temperature and/or chemical potential, one would expect that asymptotic freedom 
should make the deconfined phase of QCD accessible by analytical methods such as per- 
turbation theory JIJ. As is well known, the nonperturbative magnetostatic sector of QCD 
at high temperature produces a barrier to perturbation theory, which in the case of the 
thermal pressure occurs at order <? 6 T 4 , where g is the strong coupling constant H, [J. How- 
ever, already well below the order where this problem occurs, perturbation theory requires 
resummations of collective phenomena such as screening. While these resummations have 
a well-defined expansion parameter, the resulting perturbation series exhibit surprisingly 
poor convergence behaviour so that certain further resummations seem to be necessary to 
make perturbative results useful at interesting temperatures. This problem is in fact not 
specific to QCD, nor to gauge theories in general, but also occurs in the simplest models 
such as scalar field theory. The case of O(N) cj) 4 theory, which has a solvable (in fact ex- 
tremely simple) ./V — > oo limit has been used in Ref. 0] to analyse the properties of thermal 
perturbation series and as a toy model for recent proposals of particular partial resumma- 
tion methods such as HTL screened perturbation theory [6], [?J and approximations based 
on 2PI ^-derivable schemes ^ ^, [l(], 11]. In Ref. Jl^] 1 , the large flavour- number (Nt) limit 



of QCD has been proposed as a more interesting testing ground for various methods to 
overcome the difficulties with thermal perturbation theory, because unlike 0(N — * oo) 



1 For corrected numerical results see Ref. jl3| or the newest hep-ph version of Ref. 
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</> 4 theory large- iNfy QCD, though essentially Abelian in its remaining interaction exhibits 
many relevant phenomena such as momentum-dependent screening and damping. 2 

An exactly solvable theory which involves fermions is however also of interest with 
regard to the possibility of exploring the effects of finite chemical potential, both with 
respect to the inherent problems of thermal perturbation theory, and also beyond. The 
main nonperturbative method to investigate real QCD is certainly lattice gauge theory, 
where recently important progress has been made to also cover finite chemical potential 



1 15, [lg, [L7], [1^, 16, gy], but the extrapolation to larger chemical potential and smaller 



temperature remains uncertain. 

In this work we present the extension of the results of Ref. [12, 13] for QCD in the limit 
of large Nf to cover the entire range of temperature and chemical potential for which the 
problem that large- Nf theory only exists with a cut-off below the scale of the Landau pole 
remains negligible numerically. This is indeed the case when temperature and chemical 
potential are sufficiently below the required cut-off, which can be made exponentially large 
at small effective coupling. 3 

In particular, we obtain the exact large- Nf result for the thermal pressure and a number 
of derived quantities such as quark number susceptibilities at zero chemical potential and 
the entropy at small temperatures. We use these results to compare with known results 



from thermal perturbation theory pl| , 22, g^] obtained at small chemical potential where 
dimensional reduction ^5], |2(| is applicable. Our numerical accuracy turns out to be 
sufficiently high to permit the verification of e.g. a recent three-loop result of Vuorinen 



1 27] on quark number susceptibilities as well as a numerical coefficient in the pressure at 



zero temperature obtained long ago by Freedman and McLerran 2J, We are moreover 
able to extract a number of perturbative coefficients at order g e that are not yet known 
from analytical calculations. The comparison of the perturbative results with the exact ones 
allows us to investigate the convergence properties and ambiguities of thermal perturbation 
theory in some detail. 

For the case of both non-zero temperature and non-zero chemical potential, we are 
able to explore the range of validity of a scaling behaviour noticed recently in lattice 



calculations by Fodor, Katz, and Szabo [16| at moderately large chemical potential and 



find that it breaks down rather abruptly at fi q > irT, which points to a presumably generic 
obstruction for extrapolating data on the equation of state from small to high chemical 
potential. 

At small temperatures T <C fi, we find a region which is dominated by non- Fermi- liquid 
behaviour. There the usual linearity of the entropy in T is replaced by a nonmonotonic 
behaviour, which leads to an excess of the entropy over its free-theory value up to a certain 
value Tf fi which depends on the strength of the coupling. Thus, while large- Nf QCD does 
not have much in common with the rich phase structure of real QCD, it allows one to study 



2 A similarly complicated but purely scalar field theory in 6 dimensions, which can also be solved in 
something similar to the large- Nf has been studied in Ref. |Jw|, but being a scalar theory with cubic 
interactions it involves instabilities which render a comparison with QCD impossible. 

3 The same issue arises and has been discussed before in the exactly solvable large- iV limits of the scalar 
models of Ref. jj Q. 
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aspects of non- Fermi-liquid behaviour |3(], [32fl . which have recently turned out to be of 
relevance to the colour superconducting phases of QCD [33, 34, [3^, 36]. 



2. Recapitulation of large-A^ QCD 

In Ref. the thermal pressure of hot QCD with a large number of fermions Nf ^> 

N c ~ 1 was calculated exactly at next-to-leading order (NLO) in a large Nf expansion. In 
this limit, the gauge coupling g 2 is taken to zero such that a finite value of g 2 Nf of order 1 is 
obtained. For the two cases of massless QCD and ultrarelativistic QED the theory can then 



be solved exactly, i.e. to all orders in the effective coupling which following Ref. [12, 13] 
we define as 

9 2 N f 



9cs 



QCD. 



g 2 N f , QED. 



(2.1) 



At leading order in Nf the renormalization scale dependence is determined exactly by 
the one-loop beta function according to 



1 



1 
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This implies a Landau singularity. Following Ref. [12] we define the Landau scale Al 
such that the vacuum gauge field propagator diverges at Q 2 = A£, which leads to 



A, 



MMse 



5/6 67r 2 /g2 ff ( ms ) 



(2.3) 



The presence of a Landau singularity means that there is an irreducible ambiguity associ- 
ated with the UV completion of the theory, but in the thermal pressure this ambiguity is 
suppressed by a factor (max(T, //)/Al) 4 . 

The thermal pressure in the large- Nf limit down to order iV9 is diagrammatically given 
by an undressed fermion loop, which is of order Nj, and by a gauge boson loop with an 
arbitrary number of fermion loop insertions (plus corresponding counterterm insertions), 
which are all of order Nf because g 2 Nf ~ 0(1)- Standard Schwinger-Dyson resummation 
in the gauge boson loop includes all diagrams there are to next-to- leading order, Nf, which 
is the order we shall be interested in. 

In the imaginary-time formalism the thermodynamic pressure is obtained by perform- 
ing a sum over Matsubara frequencies, which can be replaced by a contour integral, and 
subtracting the vacuum contribution. The result for the NLO (Nf) pressure is [12] 
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Figure 1: Exact result for the large- Nf interaction pressure P — Pq normalized to N g {TT 2 T 2 + fi 

7 



as a function of g 2 s (/zms) with /ij^ s = tt 2 T 2 + fi 2 , which is the radial coordinate, and <f> = arctan — . 



where n b (uj) = l/(e w/T - 1), II 

vac is the vacuum part of the gauge-boson self energy, 

n&CQ) = (^ U Q 2 -Q»Q V ) |) , (2.5) 

and Ht and IT^ are the two independent structure functions in the thermal self energy as 
given explicitly in Ref. These cannot be given in closed form (except for their 

imaginary parts Jl3| ), but can be represented by one-dimensional integrals involving the 
fermionic distribution function nf. 

All that is needed to generalize to non- vanishing chemical potential ^ is to include the 
latter in the fermionic distribution function appearing within the self-energy expressions 
TIt and II ^ according to 

n ^ T ^ ) = l ( e (fc-d/r + i + efrf/0/r + i ) • (2 ' 6) 

When evaluating the integrals above exactly by numerical means, we can safely inte- 
grate parts proportional to n& in Minkowski space, since those are exponentially ultraviolet 
safe. For terms without n& more care is required, because the expressions are potentially 



logarithmically divergent, unless a Euclidean invariant cutoff is introduced [12]. As in 
Ref. [^, [H|] we apply a cutoff and stop the d 4 Q integration at Q 2 = aA^, varying the value 
of a between 1/4 and 1/2 to estimate the irreducible ambiguity. 

3. Results and discussion 

In Fig. |] we display our exact results 4 for the interaction pressure P — Pq oc N® , where the 



Tabulated results will be made available on-line at tittp : //hep . itp . tuwien. ac . at/~ipp/data/ 
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ideal-gas limit 
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(3.1) 



has been subtracted, for the entire fi-T plane (but reasonably below the scale Landau pole). 
For this we introduce an angle cf> = arctan — and encode the magnitudes T/Al and /u/Al 
through the running coupling g 2 s (p-Ms) with p^ s = ir 2 T 2 + fj 2 according to (2.3). 

We found that the ambiguity arising from the presence of a Landau pole reaches 
the percent level for g 2 s ft 28, where Al/ y tt 2 T 2 + fi 2 ;$ 19. At larger coupling (corre- 
sponding to larger T and/or /i), this ambiguity grows rapidly and will be shown in the 
two-dimensional plots below by a (tiny) red area. 

In the following we shall compare the exact large- Nf result with known results from 
perturbation theory at high temperature and small chemical potential, where dimensional 
reduction is an effective organizing principle, and with results at zero temperature, where 
dimensional reduction does not apply. We also investigate to what extent quark number 
susceptibilities at vanishing chemical potential determine the behaviour at larger chemical 
potential. 

3.1 Pressure at small chemical potential 

The perturbative result for the thermal pressure of hot gauge theories with fermions has 
been obtained to order g° at zero chemical potential in Ref. [^, 22 1 and, using effective 
dimensionally reduced field theory, in Ref. |37| ]. 

In the large- Nf limit, dimensional reduction in fact gives a free theory at order N®. 
Below order g 7 , its pressure contribution is simply Pqr = N g m^T '/ '(127T 2 ), where ttle 
is the Debye mass obtained by integrating out hard momentum modes to the required 
perturbative order. (At and beyond order g 7 , the effective theory requires higher-derivative 
kinetic terms.) 

Dimensional reduction continues to work for small chemical potential fj, < T. The 
result to order g^ s for QCD in the large- Nf limit as available in the literature reads 
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where the terms oc g 4 ff and involving [i have recently been computed by Vuorinen |27], . 
The contribution to order ff 3 follows from the NLO correction to the effective-field-theory 



parameter m E , computed at finite \i in Ref. [23 
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with the function 



;°° dp ( 1 - eP 1 1 

-oo 



-oo P V 1 + eP eP+nx + 1 e~P +nx + 1 1 
= + = -2 7 -41n2-2Re^(^(l + «x)), (3.4) 

where VK- 2 ) is the digamma function and = ip(z) + is notation introduced by 

Ref. |38|]. 5 For small x this function can be expanded as 

V(x) = 4 E(-l) n (l " ^+i) C(2n + 1) x 2 " , (3.5) 

evidently with a radius of convergence of 1, which corresponds to /i = ttT. 

Because the dimensionally reduced theory is free, the only nonanalytic terms in g^ s 
stem from the "plasmon term" oc m%. In particular there are no logarithms ln(g) which in 
finite-Ay QCD appear at and beyond order g 4 and which have recently been determined 
even to order g 6 in Ref. |M| . 



The numerical result for [i = 0, obtained before in Ref. [12, 13], is shown in Fig. but 
now as a function of 5^ (/2ms = ttT). Also given are the respective values of log 10 (AL/7rT). 

For small coupling the agreement of our numerical results with perturbation theory is 
sufficiently accurate that it permits a numerical extraction of the coefficients to order g^ s 
which are not yet known analytically. 

Eq. ( |2.2j ) dictates that the <7g ff -term in the pressure at \i = has the form 

1 



= ^Vr 4 [c fi + ioin 2 ^ 



g%,l*=0 V 4tt / 



irT 



( 1 + _ 4641n2 _ ^ + 16 0^\ [D ~ mS 



81 V 5 C(-3) C(-l) / 



(3.6) 



and by least-square fitting we obtain numerically the estimate C% = +20(2). 

In real, finite-iVj- QCD this result corresponds to the coefficient involving Nj in the 
g 6 T 4 term of the pressure, which is nonperturbative in its purely gluonic contribution 



In Ref. 0, |l| the convergence of successive perturbative approximations to order g^ s 



has been studied, with the result that there are large renormalization scale dependences 
beyond g^ s ~ 4. Fixing this scale dependence by the requirement of "fastest apparent 
convergence" (FAC) in the m 2 E parameter of dimensional reduction leads to /2ms = Afac = 
7re 1 ' /2_7 T. This choice leads to fairly accurate results up to g 2 ff ~ 9. 

At still larger coupling the exact result for the pressure at zero chemical potential has 
the remarkable feature of bottoming out and tending toward the ideal-gas limit. Before 



5 The closed-form result for T> can be obtained from Ref. ]38|] after identifying 

V{x) = (4tt) 2 [j 2 ( M =-7rTx-) -22(/x = 0)] 

(cf. Eq. (5.4) of [Hi and the definition (A. 2) of Q) and using the result (B.71) for Ti along with the 
definitions (2.15), (2.16), and (2.20) of H. 
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Figure 2: Exact result for the interaction pressure at zero chemical potential as in Ref. |12|, |l3fl 
but as a function of g^. s (fiMS = kT) or, alternatively, log 10 (AL/7rT). The purple dashed line is the 
perturbative result when the latter is evaluated with renormalization scale /2ms = /2fac = Tre 1 / 2 r T; 
the blue dash-dotted lines include the numerically determined coefficient to order g® s (with its 
estimated error) at the same renormalization scale. The result marked "g^ s = g® s " corresponds 
to choosing /2ms such that the order-g® s coefficient vanishes and retaining all higher-order terms 
contained in the plasmon term oc m^. In this and the following plots the (tiny) red band appearing 
around the exact result at large coupling displays the effect of varying the cut-off from 50% to 70% 
of the Landau scale Al- 



it begins to exceed the latter, however, the Landau pole starts to influence the result 
noticeably. This is displayed in the figures by the tiny red area around the "exact" result, 
which represents the effect of varying the cut-off from 70% to 50% of the Landau scale 



Al- (See Ref. [13| for the effect of varying the cut-off independently in the Minkowski and 
Euclidean parts of the numerical integrations.) 

Including our numerical estimate of the (^-coefficient and using /2fac further improves 
the perturbative result so that it remains accurate up to g% s ~ 16. The agreement with 
the exact result can even be further improved by fixing the renormalization point such 
that the g® s coefficient vanishes and keeping all orders of the odd terms in g e s by leaving 
the plasmon term oc m| unexpanded in g^ s . The result of this procedure is indicated 
by the gray area in Fig. |2[ The apparent success is in line with the recent observation 



in Ref. [4C] that keeping the parameters of the dimensionally reduced theory unexpanded 



greatly improves the convergence of thermal perturbation theory. 

3.2 Quark number susceptibilities 

3.2.1 Linear quark number susceptibility 

The (linear) quark number susceptibility is defined as the first derivate of the quark number 
density M with respect to chemical potential, 

dM d 2 P 

x= w = w (3 - 7) 
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Figure 3: The interaction part of the quark number susceptibility at /i = compared with strict 
perturbation theory to order and g® s , respectively, with renormalization scale varied about irT 
by a factor of 2. 



Fig. H displays the exact large- Nf result for the interaction part of % at zero chemical 
potential as a function of g c s (or alternatively log 10 (AL/vrT)). Similar to the thermal 
pressure, the result is nonmonotonic, but the minimum already occurs at g 2 s (irT) ~ 8.6, 
and the free-theory value is recovered at g 2 s (iTT) ~ 22.5, where the Landau ambiguity is 
still well under control since Al/T ~ 100 at that coupling. 

The perturbative (dimensional reduction) result can be read from the linear term in 



\x of (|3.2|) and gives 
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(3.8) 



The coefficient of has only recently been obtained in [27| in a three- loop calculation. 
We can confirm its closed- form value by a numerical fit, which gives agreement with an 
accuracy of 2 x 10 -4 , thus providing a good check on both our numerics and the analytical 
calculations of J27]]. This level of accuracy allows us to also extract the order-g 6 term as 
(for /i M s = ttT) 



NgT 2 



-4.55(9) x 



(9cS\ 
V4vr/ 



(3.9) 



In Fig. H| we show the perturbative results to order g^ and g^ s , varying the renor- 
malization scale about ttT by a factor of 2 (now without the improvement of keeping 
effective-theory parameters unexpanded). The value /2ms = is in fact close to JIyac 
where it makes no difference whether m 2 E is kept unexpanded or not. We find that the 
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Figure 4: The interaction part of the higher-order quark number susceptibility d^(x— Xo) at /i = 
compared with strict perturbation theory to order g^ s and g® s , respectively, with renormalization 
scale varied about irT by a factor of 2. The coloured bands of the (^-results cover the estimated 
error of the numerically extracted perturbative coefficients. 



quality of the perturbative result for the susceptibility is comparable to that observed in 
the pressure, with /2ms = being close to the optimal choice. 



3.2.2 Higher-order quark susceptibility 

We have also computed explicitly the higher-order susceptibility 9 4 P/9^ 4 | M =o (which has 



recently been investigated in lattice QCD with Nf = 2 in Ref. |2( 

Our exact result in the large- Nf limit is shown in Fig. |I[ In this quantity, we find that 
the nonmonotonic behaviour observed above in the pressure and the linear susceptibility 
is much more pronounced. The minimum now occurs at g 2 s ~ 3.7, where perturbation 
theory is still in good shape, and the free-theory value is exceeded for g 2 s ^ 9. Using (| 
we find to order g 5 



X ~ Xo 



dfi 2 Ng 
35, 



<9 4 P-P, 



M =o d/x 4 Ng 
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+ 0( 5e 6 ff) 



(3.10) 



In a first version of this paper we have extracted the coefficient appearing at order 
g 4 ff numerically as C4 = —7.02(3). In the meantime, the complete ji dependence of the 



dimensional reduction result to order g 4 ff has been worked out in Ref. [38] from where one 
can obtain the exact result 



C 4 = -— - 12 In 2 + -C(3) 
12 6 SV ; 



-6.9986997796998 . 



(3.11) 
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Figure 5: Deviation from the scaling observed in Ref. |16[ in lattice QCD for small chemical 
potential in the quantity SP = P(T, n) ~ P(T, 0) - ^x\^o(p 2 + A* 4 I '(2ir 2 T 2 )) (full lines) and in 5P = 
P(T,n) - P(T,0) - £ X \n=o - 70l=o (dashed lines), both normalized to pgfo = N g n 2 T 4 /45, 
for 5e 2 ff (7rT) = 1,4,9,16. 



The complete agreement with Ref. |^] provides on the one hand an independent check 



on the correctness of the 3- loop calculations of 1 38] and on the other hand a check on the 
accuracy of our numerical analysis. 



Using ( |3.11 ) we can extract the term of order g® s in ( p. 10 ) as —39(l)g® s /(l28ir 8 ) for 



fi = ttT. The perturbative results to order <Aj and to order g^ s are compared with the exact 
result in Fig. |j. This shows that the accuracy of the perturbative result again improves 
by going from order g^ s to order g® s , but the renormalization scale dependence increases 
sharply at large coupling. 6 

3.2.3 Pressure at larger chemical potential from susceptibilities 

With regard to the recent attempts to explore QCD at finite chemical potential by means 



of lattice gauge theory []16j, |17], |18|, |19f| , it is of interest how well the pressure at larger 
chemical potential can be approximated by the first few terms of a Taylor series in \j? . 



In Ref. [16] it has been observed that the ratio of AP = P(T,fi) — P(T,fi = 0) over 
the corresponding free-theory quantity APo is practically independent of \jl for the range 
of chemical potentials explored. This is also realized when quasi-particle models are used 
for a phenomenological extrapolation of lattice data |^l|, Q in a method introduced by 
Peshier et al. p3|]. 

In Fig. [| we show the deviation from this "scaling" at higher values of /u/T by con- 
sidering the quantity 5P = P(T,/i) - P(T,0) - \x\^{p 2 + ^ 4 /(2vr 2 T 2 )) divided by 
pfrf* = jV 9 7r 2 T 4 /45. The combination + n 4 /(2ir 2 T 2 )) appearing therein is such that 
a replacement of P and \ by their interaction-free values Pq and xo makes 5P vanish 



6 The size of the scale dependence can in fact be reduced somewhat by keeping odd powers of itie without 
expanding out the correction in m%. 
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identically. (As can be seen from the above perturbative results, 5P also vanishes for the 
leading-order interaction parts oc <7g ff .) In the exact large- Nj results of Fig. || we observe 
that for coupling g^ s < 4 the deviation 5P is at most a few percent of -PjsjLO ^ or /V^ ~ ^ 
but it rapidly grows for fx/T > ir. This is in fact also nicely illustrated by the 3-dimensional 
plot of the pressure in Fig. |l|, which has a rather conspicuous kink at <j) = 45° corresponding 
to = ttT. 

It turns out that including the exact higher-order susceptibility at /j, = does not lead 
to a better approximation of the pressure at larger chemical potential. The dashed lines 
in Fig. | correspond to 5P = P{T, n) - P(T, 0) - *fxU=o - ir0U=o- While this slightly 
improves matters at small /i/T, it results into even quicker deviations for larger /i/T. 

It is of course impossible to say whether this behaviour would also appear in real 
QCD, but since it occurs already at comparatively small g e fi in the large-iVy limit, where 
the peculiar nonmonotonic behaviour of the pressure as a function of g e d does not yet arise 
(the minimum in the normalized interaction pressure occurs at g^ s (irT) ~ 14), 7 it may be 
taken as an indication that extrapolations of lattice data on the equation of state from 
small chemical potential to large fj,/T are generally problematic. If anything, real QCD 
should be more complicated because of the existence of phase transitions which are absent 
at NLO in the large- iVj limit. 

3.3 Pressure at zero temperature 

Our exact result for the thermal pressure at zero temperature and finite chemical potential 
is given in Fig. || as a function of g^ s (/2ms = A 4 ) ■ I n contrast to the pressure at zero chemical 
potential and finite temperature, the interaction pressure divided by /Z 4 is monotonically 
decreasing essentially all the way up to the point where the Landau ambiguity becomes 
noticeable. 

The thermal pressure at zero temperature and large chemical potential for QED and 



QCD has been obtained to order g A long ago by Freedman and McLerran |28| , 29]. At this 
order, there is a non-analytic zero-temperature plasmon term oc g ln(g), whose prefactor is 
known exactly, but the constant under the logarithm only numerically. The transposition 
of their result, which has been obtained in a particular momentum-subtraction scheme, to 
the gauge-independent MS scheme can be found in Refs. |45| ]. The large- Nj limit of this 
result reads 



P-P 



2 



T=o 32ir 4 



1„M _l ln tMS_c 4 
2tt 2 3 n 



+0(g e cfi lng cS ) (3.12) 



1287T 



and involves one of the numerical constants computed in Ref. 



~ 79 7T 2 7 log 2) 2b , 

C 4 = — w 0.536, (3.13) 

18 3 3 3 v ' 



where b has an integral representation, given in Eq. (II. 3. 25) of Ref. |28j, that apparently 



cannot be evaluated in closed form. Fixing an obvious typo in Eq. (II. 3. 25) of [28], b 
7 Peshier has recently argued that the strong-coupling behaviour of large Nj QCD will be relevant 

for real QCD at most in the coupling range where the normalized pressure decreases with g c a . 

8 Comparison with Eq. (II. 3. 24) shows that there is a missing exponent 2 after the second set of large 

round parenthesis in Eq. (II. 3. 25) of j28|. 
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Figure 6: The interaction part of the pressure at zero temperature and finite chemical potential 
as a function of g^ s (fiMS — A 4 ) or, alternatively, log 10 (AL/Ai), compared with the perturbative result 
of Freedman and McLerran p8j, |2£j to order g* s , and our numerically extracted order- g® s result, 
both with renormalization scale in the perturbative results varied around /2ms — A 1 by a factor of 
2. The coloured bands of the (^-results cover the error of the numerically extracted perturbative 
coefficients. 



can however be easily evaluated numerically to higher accuracy than that given in |28J as 
b = -1.581231511 . . ., which leads to C 4 = 0.5358316747 . . . . 

The accuracy of our numerical results is sufficiently high to confirm the correctness of 
the result for C4 with an accuracy of ~ 2 x 10~ 4 . With the knowledge of the exact value 
of C4 we can also extract, with lower precision, the next coefficients at order g® s , which 
again involve a logarithmic term: 

P-Po 



T=0 



9 2 eS 

32vr 4 



In 



9 2 es 



2vr 2 
3.18(5) - 



2 . /Ums 
— in 

3 /i 

3 \x 



C 4 



1287T 6 



In 



9cs 



16 2 /XMS 

9 n 



ln 



2vr 2 



(3.14) 



3.4(3) 



9eff 



2048vr 8 



+ 



In Fig. ^ we also study the renormalization scale dependence and apparent convergence 
of the perturbative result. We have varied /2ms about \i by a factor of 2, and it emerges 
that the larger values are somewhat favoured. 

At low temperature T <C [i, dimensional reduction does not occur. If one nevertheless 
considers the effective- field-theory parameter m 2 E of ( |3.3|) in this limit, one finds that the 
function T>(x) therein approaches —2 (ln 2x + 7), so that the T 
reads 



limit of m? E exists and 



ft 



7T 2 



6vr 2 



ln 



/^MS 

2fi 



+ 0{g\ 



eg) 



(3.15) 



Fastest apparent convergence applied to this quantity would suggest /iMS = 2e2/i 3.3/i. 
This turns out to be not as good as the choice of 2/i, though slightly better than /2ms = A 4 - 
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3.4 Entropy at small temperatures and non-Fermi-liquid behaviour 

The effect of small temperature on the pressure at nonzero chemical potential can be 
studied in terms of the entropy 

(3-16) 



vM7„ 

from which the specific heat can be derived. (Both, the entropy and the specific heat vanish 
in the zero-temperature limit, and the various kinds of specific heat have the same small-T 
behaviour up to terms ~ T 2 .) 

At small temperature T <C \i one might expect the contributions involving the Bose- 
Einstein distribution in the thermodynamic potential ( |2.4| ) to be negligible compared 
to the "non-n?," contributions. 

Considering the latter contributions first, we find that the corresponding part of the 
entropy vanishes linearly 

<S Q on-n b ~> T(T for T -> (3.17) 

with 

^ + 0(gt {i lng efi ). (3.18) 



G — °"0 _ d 2 P aon -n b — Po 



N gf i 2 dT 2 N g n 2 



T=o 8ir 2 



The coefficient at order g 2 s is in accordance with the g 2 s part of the strictly perturbative 
result for the pressure, which is also known as the exchange term ffl and which coincides 
with the g 2 s part of (fg). 

At small coupling, we can also extract the order-g 4 ln(<?) corrections to a from S non . nb 
numerically with the result 



1 «\ & 



Ngfi 2 ' 5 off 32vr 4 



- In ^ - 0.328(1) x In M + 0.462(5) 
3 \i 2ir 2 



+ ... (3.19) 



The exact result for S non _ nb is given by the dash-dotted lines in Fig. [7| for g 2 s (flMS = 
fi) = 1, 4, and 9 and < T/fj, < 0.15. In this range of temperatures, S non . nb is well 
approximated by the linear term ( |3.17j ). 

Numerically evaluating also S nb , i.e. the contributions to the entropy following from 
the parts of fl2.4|) which involve n^, we find that these cannot be neglected at small temper- 
atures. As shown in Fig. 0, for sufficiently small T j ' ji, S nb is positive and even dominates 
so that the total result for the entropy turns out to exceed its free-theory value for a certain 
range of T/fi, which gets larger with increasing g 2 s . 

The largest part of the positive and nonlinear contributions at small T / fi in fact comes 
from the transverse vector-boson modes. At small frequencies, these are only weakly dy- 
namically screened (and completely unscreened in the static limit because of gauge invari- 
ance). As has been discussed in a variety of contexts (nonrelativistic and relativistic QED 
as well as colour superconducting QCD) in Refs. |3(], [H], [32|, [33|, |34], |35|, 36 1, this fact gives 
rise to non-Fermi-liquid behaviour at sufficiently small temperature. 

A particular consequence is the appearance of anomalous contributions to the entropy 



as well as specific heat. In Refs. |30|, 31, 33] these contributions have been found to be of the 



order — g 2 s [i T In T. This has however been questioned recently by the authors of Ref. p6| 
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Figure 7: The interaction part of the entropy at small T//i for g^ s (p,yis = /-*) = 1, 4, and 9. 
The "non-nb" contributions (dash-dotted lines) are negative and approximately linear in T with 
a coefficient agreeing with the exchange term cx g^ s in the pressure at small coupling; the "rib" 
contributions (dashed lines), which are dominated by transverse gauge boson modes, are positive 
and nonlinear in T such that the total entropy exceeds the free-theory value at sufficiently small 
T/fi. 



who found a different behaviour which in the large- Nf limit would imply — g 2 ff T 3 mT 
as the dominant non- Fermi-liquid contribution, with TlnT contributions appearing only 
at higher order in the coupling. Because in both of these apparently conflicting results 
there are also non-anomalous and not yet determined terms such as g^ s fj?T, we cannot 
discriminate with certainty between the two. On the other hand, we could verify that the 
coefficient of the —TlnT contribution as given in 32] has about the correct magnitude 
to permit a reasonable fit to our exact results, favouring in fact a leading power of g^ s in 
this contribution for T j fi < <^ ff /(47r 2 ) (where the analysis of Ref. [36| may well cease to be 
applicable) . 

We intend to investigate this matter however in more detail in a separate work. 



4. Conclusion 

We have extended the previously obtained exact result for the pressure of hot QCD in 
the limit of large flavour number to finite chemical potential. For small coupling we have 
been able to confirm numerically a number of previously calculated perturbative coeffi- 
cients with great accuracy, and have even obtained numerical values for a few perturbative 
coefficients at order g G which were not yet known. At larger coupling we have studied 
the applicability of the perturbative results, their renormalization scale dependence and 
apparent convergence at the highest available order in the coupling. 

For the region of simultaneously non-zero temperature and chemical potential, we have 
observed a comparatively weak dependence of the pressure on the chemical potential for 
[i < ttT, which is correspondingly well described by the quark number susceptibility at 
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zero chemical potential. At fx ~ ttT, we instead observed a rather abrupt change leading 
to a breakdown of the small-/i scaling. We suspect that also in real QCD, where a simple 
scaling determined by the quark number susceptibility has been observed recently in lattice 
QCD [16 1, this may be similarly limited. 

At small temperatures T< /jwe found non- Fermi-liquid contributions to the entropy. 
As a result, for a certain range of temperature, entropy and specific heat show strong 
deviations from linearity in T and also an excess over their free-theory values. The exactly 
solvable large-TVy limit of QED and QCD evidently allows one to study the phenomenon 
of non-Fermi-liquid behaviour beyond perturbation theory. A more detailed investigation 
will be the subject of a future publication. 
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